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Abstract 



We investigate three models for the charge symmetry breaking (CSB) of 
the nucleon-nucleon {NN) interaction (based upon p-uj mixing, nucleon mass 
splitting, and phenomenology) that all reproduce the empirical value for the 
CSB of the ^Sq scattering length {Aacss) accurately. We reveal that these 
models make very different predictions for CSB in ^Pj waves and examine 
the impact of this on some observable quantities of A > 3 nuclear systems. 
It turns out that the ^H-'^He binding energy difference is essentially ruled by 
AacsB and not very sensitive to CSB from P waves. However, the Coulomb 
displacement energies (which are the subject of the Nolen-Schiffer anomaly) 
receive about 50% of their CSB contribution from A^A^ partial waves beyond 
^5*0. Consequently, the predictions by the various CSB models differ here 
substantially (10-20%). Unfortunately, the evaluation of the leading Coulomb 
contributions carry a large uncertainty such that no discrimination between 
the competing CSB models can presently be made. To decide the issue we 
suggest to look into nuclear few-body reactions that are sensitive to CSB of 
the nuclear force. 
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I. INTRODUCTION 



By definition, charge independence is invariance under any rotation in isospin space. 
A violation of tliis symmetry is referred to as charge dependence or charge independence 
breaking (CIB). Charge symmetry is invariance under a rotation by 180° about the y-axis in 
isospin space if the positive z-direction is associated with the positive charge. The violation 
of this symmetry is known as charge symmetry breaking (CSB). Obviously, CSB is a special 
case of charge dependence. 

CIB of the strong A^A^ interaction means that, in the isospin T = 1 state, the proton- 
proton (Tz = +1), neutron-proton (Tj = 0), or neutron-neutron {T^ = —1) interactions are 
(slightly) different, after electromagnetic effects have been removed. CSB of the A^A^ interac- 
tion refers to a difference between proton-proton (pp) and neutron-neutron (nn) interactions, 
only. 

Charge asymmetry is seen most clearly in the ^5*0 scattering length. The latest em- 
pirical values of the neutron-neutron [nn) singlet effective range parameters (corrected for 
electromagnetic effects) are 

al = -18.9±0.4fm[l,|, r^„ = 2.75 ± 0.11 fm §. (1.1) 

This should be compared to the corresponding proton-proton (pp) values |0: 

= -17.3 ± 0.4 fm, = 2.85 ± 0.04 fm. (1.2) 

The implication is that the singlet effective range parameters break charge-symmetry by the 
following amounts, 

AacsB = <p - C = 1.6 ± 0.6 fm, (1.3) 
ArcsB = r^p- = 0.10 ± 0.12 fm. (1.4) 

The current understanding is that — on a fundamental level — the charge dependence of 
nuclear forces is due to a difference between the up and down quark masses and electromag- 
netic interactions among the quarks. As a consequence of this — on the hadronic level — major 
causes of CIB are mass differences between hadrons of the same isospin multiplet, meson 
mixing, and irreducible meson-photon exchanges. For recent reviews on charge dependence, 
see Refs. §-§. 

Neutral mesons with the same spin and parity, but different isospin, may mix due to the 
up-down quark mass difference and electromagnetic interactions; that is, the Hamiltonian 
responsible for the mixing, Hm, has a strong and an electromagnetic part, 

Hm = Hstr + Hem ■ (1-5) 

The most prominent case is p^-u mixing that is observed in the annihilation process e^e^ — *■ 
TT+TT^ from which the mixing matrix element has been extracted to be 

{p^\Hm\uj) = -0.00452 ± 0.0006 GeV^. (1.6) 

The charge-asymmetric nuclear force created by this process is displayed in Fig. ^ Coon 
and Barrett showed that this mechanism alone can explain the entire empirical CSB of 
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the singlet scattering length, Eq. (|1.3| ). Other examples of meson mixing are tt-t] and ir-r]' 
mixing which, however, generate only negligible effects [0. 

In recent years, the process of Fig. |l]has been subjected to some criticism. Note that the 
matrix element, Eq. (|1.6| ), is extracted on-shell, i. e., for k"^ = m^, where k denotes the four- 
momentum of the p meson and rrip the p mass. However, in the A^A^ interaction. Fig. |1], the 
relevant k^ are space-like (less than zero). Using a quark model for Hgtr, Goldman, Hender- 
son, and Thomas find a substantial fc^ dependence which is such that the contributions 
of Fig. in would nearly vanish. Similar results were reported in subsequent papers P-[TT[]. On 



the other hand. Miller ^ and Coon and coworkers have advanced counter-arguments 
that would restore the traditional role of p-u exchange. The issue is unresolved. Informative 
summaries of the controversial points of view can be found in Refs. [p|,p!^,p!^ . 

We now turn to another basic source for CSB of the nuclear force, namely, nucleon mass 
splitting. The most trivial consequence of nucleon mass splitting is a difference in the kinetic 
energies: for the heavier neutrons, the kinetic energy is smaller than for protons. This raises 
the magnitude of the nn scattering length by 0.25 fm as compared to pp. Besides this, 
nucleon mass splitting has an impact on all meson-exchange diagrams that contribute to 
the nuclear force. For example, there are the one-boson exchange (OBE) diagrams. Fig. |^, 
which are affected by only a negligible amount. However, a sophisticated and realistic 
meson model for the nuclear force |15| goes beyond single meson exchange and includes 
irreducible diagrams of two-boson exchange which generate substantial CSB from nucleon 
mass splitting. The major part of the CSB effect comes from diagrams of 27r exchange with 
A^A intermediate states. Fig. |^. In fact, these diagrams can fully explain the empirical CSB 
splitting of the singlet scattering length |l6| , p!7| . 

Finally, for reasons of completeness, we mention that irreducible diagrams of vr and 7 
exchange between two nucleons create a charge-dependent nuclear force. Recently, these 
contributions have been calculated to leading order in chiral perturbation theory |jl8|. It 
turns out that to this order the 717 force is charge-symmetric (but does break charge inde- 
pendence). 

The bottom line then is that we have two CSB mechanisms in hand, each powerful 
enough to fully explain the charge asymmetry seen in the singlet scattering length. This 
state of affairs challenges the question: Which of the two mechanisms is nature really using? 
To answer the question, the ^5*0 state at zero energy is obviously of no use since this can be 
described equally well by both mechanisms. Thus, the answer, if any, can only come from 
energies Tiat > and/or states with L > 0. It is quite possible that the predictions by the 
two mechanisms are very different in states other than ^Sq. The 27r exchange contribution 
to the nuclear force (where the CSB effect due to nucleon mass splitting mainly comes from) 
is chiefly a central force of intermediate range, while vector meson exchange (involved in p-Lj 
mixing) creates a large spin-orbit component and a short-ranged central force. This may 
create large differences in triplet P waves where the signature of spin-orbit forces is most 
pronounced. If large differences between the models occur for L > 0, then each mechanism 
may have a characteristic signature on observables that are sensitive to CSB from P and 
higher partial waves. In Ref. it was found that about 50% of the CSB contribution 



to the Nolen-Schiffer (NS) anomaly comes from the L > two-nucleon states. Thus, 
there is the possibility that the differences between competing CSB models may show up 
as differences in the predictions for the NS anomaly. Other observables may be considered, 
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too. 

It is the purpose of this paper to look into this issue. In particular, we want to find out if 
it is possible to discriminate between the two models for CSB based upon their predictions 
for observable quantities, other than the singlet scattering length. 

Besides the two CSB mechanisms discussed, we will include in our study also a phe- 
nomenological model for CSB, namely, the Argonne Vis potential ||2l|. The three models 



will be introduced in Sec. II and applied to the ^H-^He binding energy difference and the 
NS anomaly in Sec. III. Sec. IV concludes the paper. 



II. MODELS FOR CHARGE SYMMETRY BREAKING OF THE 
TWO-NUCLEON INTERACTION 



A. p-uj mixing 

Here we will evaluate the charge asymmetric nuclear force generated by p^-u) mixing 
shown in Fig. |l| |2[. 

The coupling of p and u) mesons to nucleons is described by the following Lagrangians: 



c 



pNN 



ip) 



4M 



(2.1) 
(2.2) 



where ip denotes nucleon fields, ip meson fields, and Xj 



1,2,3) are the usual Pauli 



matrices describing isospin specifically, Ts\p) = +\p) and Tain) = —\n) with \p) a proton 
state and \n) a neutron state. Mp is the proton mass which is used as scaling mass in the 
pNN Lagrangian to make fp dimensionless. To avoid the creation of unmotivated charge 
dependence, the scaling mass Mp is used in pp as well as nn scattering. 

The first Feynman diagram of Fig. |l] leads to the following p-u potential. 



(q'A;A^|l^;^)|qAiA2) 



1 M M 
(27r)3 



3 \l^ {p''\Hjuj) gp 



xn(q, A'l) Tg 



7. + J^^^Al' - lY 



u{q, Ai 



XM(-q',A2) 7" M(-q, A2) 
[(q' - q)2 + mj] [(q' - q)2 + ml] 



(2.3) 



where M denotes the relevant nucleon mass (i. e., M = Mp in pp scattering and M = Mn 
in nn scattering), nia are meson masses, and Hp = fp/gp- We are working here in the two- 
nucleon cm. frame and use the helicity formalism. The helicity Aj of nucleon i is defined 
as the eigenvalue of the helicity operator ^cTj ■ Pi/|pi| which is ±|. In-coming nucleon 
1 carries helicity Ai and momentum q and in-coming nucleon 2 carries helicity A2 and 
momentum — q; the out-going nucleons have A'^, q' and A2, — q'. The above 'quasi-potential' 
is defined in the framework of the relativistic, three-dimensional Blankenbecler-Sugar (BbS) 
equation which is a reduced version of the four-dimensional relativistic Bethe-Salpeter 
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equation . In the BbS formalism, the four-momentum transfer between the two nucleons 
is (g' — g)^ = (0, q' — q). The square-root factor M/\/E'E (with E = + q^ and 

E' = \l M'^ + q'^) makes it possible to cast the BbS equation into a form that is identical to 
the conventional Lippmann-Schwinger equation. The Dirac spinors in helicity representation 
are given by 

«(-1.A.)^/^(^)|A.>, (2.5) 

with normalization 



u{q,X)u{q,X) = l, (2.6) 
where u = u^7°. The above amplitude, includes form factors 

with rua the mass of the meson involved and Aq the so-called cutoff mass. For more details 
concerning the formalism and the explicit evaluation of quasi-potentials of this kind, see 
Appendix E of Ref. 0. 

The second Feynman diagram of Fig. ^ yields the same as the first one and, so, the entire 
p-u potentials is 

Vp. = 2V,« (2.8) 

with V^i^^ given by Eq. (|2.3| ). Since {p^\Hm\uj) is negative, V^^ is repulsive for pp scattering 
and attractive for nn scattering with the magnitude of V^^ essentially the same in both 
cases. 

When constructing CSB A^A^ potentials, one starts with the pp potential since there are 
many pp data (and pp phase shift analyses) of high quality to constrain the pp potential. 
In this work, we use the CD-Bonn pp potential which reproduces the world pp data 
below 350 MeV lab. energy available in the year of 2000 with the perfect x^/datum of 1.01. 
The nn potential is then fabricated by adding to the pp potential a difference-potential that 
contains the entire difference between nn and pp due to p-u exchange, which is 

AVp^ = V,^{nn) - Vp^PP) ~ -W^^^\M^) (2.9) 

where V^^^ is V^^^ [Eq. ( p.3| )] with the operator replaced by 1. Note that — strictly 
speaking — Vp^{nn) is to be evaluated with M = M„ and Vp^{pp) with M = Mp. How- 
ever, if we wish to subsume both terms into one, then we have to use the same mass for 
both for which we choose M = Mn- We have tested this approximation and found that it 
affects the singlet scattering length by 10~^ fm. 

To obtain a convenient expression for AVp^, we make use of the identity 
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[(q' - q)2 + ml] [(q' - q)^ + ml 



mt 



(q' - q)2 



m. 



1 

(q' - q)2 + ml 



which allows us to write AV^^ in terms of the difference of two expressions each 
resembles single meson exchange; namely, 



, (2.10) 
of which 

(2.11) 



with 



(q'A;A'2|AKL^)|qAiA2) 



/M M {p^\Hra\uj) 



(2vr) 



E' \ E m2 



9p 9u 



XM(q', A'l) 



2M, 



M(q, Ai 



xn(-q',A^) 7^ M(-q, A2) 
., ^p[(q^-q)^]^.[(q^-q)^] 

(q' ~ q)^ + 



and 



(q'A;A^|AK(:r)|qAiA2) 



Oif [m {p^\Hm\u}) 



(27r)3 
XM(q', A'l 



E' \ E m^, — mi 



9p 9u 



K 



7p + j^vw - <iy 



M(q, Ai) 



X«(-q',A'2)7^M(-q, A2) 
^p[(q'-q)^]-^.[(q'-q)^] 



X 



(q' - q)^ + f^l 



For the masses involved, we use [E6 



Mp = 938.27231 MeV, 
M = Mn = 939.56563 MeV, 



(2.12) 



769.9 MeV, 
781.94 MeV, 



(2.13) 



(2.14) 
(2.15) 
(2.16) 
(2.17) 



and Ap = A, 



1.4 GeV. 

We choose for the meson-nucleon coupling constants. 



^ = 0.84, 



47r 



Arc 



6.1 



^ = 10, 



(2.18) 
(2.19) 

(2.20) 



and for the mixing matrix element. 
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{p^\Hra\oj) = -0.00402 GeV2, (2.21) 

to obtain 

AacsB = 1-508 fm. (2.22) 



The above mixing matrix element is consistent with the empirical value, Eq. ( |1.6|) , and the 
p parameters are identical to the ones used in the Bonn Full Model [|T5|. Concerning the 
u, the Bonn model uses g'^/4:7c = 20 which, however, would generate too much CSB when 
applied in the above p-u potential; therefore our choice Eq. ( |2.20| ). This choice could be 
justified with the argument that part of the u contribution in meson-theoretic potentials 
may be just a parametrization of short-ranged repulsion that is actually due to quark-gluon 
exchange. 

To check our calculations, we have made a comparison with the results by Coon and 
Barrett Note that these authors use very different vector-meson coupling constants as 
compared to ours. In terms of our convention for the coupling constants, Eqs. (p.l|) and 
(^ 1^, Coon and Barrett use gl/Air = 0.6, Kp = 3.7, gl/An = 5.25, and = -0.12. 



There are also other differences, like, we use the the full relativistic Feynman amplitudes 
for p^-uj exchange, Eq. ( |2.3| ), while in Ref. |Q the nonrelativistic approximation is applied. 
Moreover, Coon and Barrett use the Reid potential ||28| as their basic pp potential whereas 



we use the pp CD-Bonn Taking all these differences into account, we were able to show 
that our results are consistent with the findings of Coon and Barrett. 

The bottom line is that due the uncertainties in the model parameters, there is latitude 
of a factor of two or so in the strength of the p-uj potential. Within that latitude, it is easy 
to fit the full amount of CSB of the singlet scattering length. 



B. Nucleon mass splitting 

The difference between the masses of neutron and proton represents a basic cause for 
CSB of the nuclear force. This source of CSB effects has been explored in great detail in 



|[r?| . The investigation is based upon the Bonn model for the NN interaction [15]. Let 
us briefly summarize the results. For this we divide the total number of meson exchange 
diagrams that is involved in the nuclear force into several classes. Below, we report the 
results for each class separately. 

1. One-boson-exchange (OBE, Fig. ^ contributions mediated by 7r°(135), p°(770), 



u;(782), ao/5(980), and cr'(550). In the Bonn model [|T5|, the a' describes only the 
correlated 27r exchange in tttc — S-wave (and not the uncorrelated 2n exchange since 
the latter is calculated explicitly, cf. Fig. Charge- symmetry is broken by the fact 
that for pp scattering the proton mass is used in the Dirac spinors representing the four 
external legs [Fig. 0(a)], while for nn scattering the neutron mass is applied [Fig. @(b)]. 
The CSB effect from the OBE diagrams is very small. 

27r-exchange diagrams. This class consists of three groups; namely the diagrams 
with NN, NA and AA intermediate states, where A refers to the baryon with spin 
and isospin | and mass 1232 MeV. The most important group is the one with A^A 
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intermediate states which we show in Fig. ^. Part (a) of Fig. 2 apphes to pp scattering, 
while part (b) refers to nn scattering. When charged-pion exchange is involved, the 
intermediate-state nucleon differs from that of the external legs. This is one of the 
sources for CSB from this group of diagrams. The 27r class of diagrams causes the 
largest CSB effect. 

3. 7rp-exchanges. Graphically, the irp diagrams can be obtained by replacing in each 27r 
diagram (e. g., in Fig. |^) one pion by a p-meson of the same charge state. The effect 
is typically opposite to the one from 27r exchange. 

4. Further Svr and An contributions {ira + TTUj). The Bonn potential also includes some 
Svr-exchanges that can be approximated in terms of na diagrams and 47r-exchanges of 
Tiuj type. The sum of the two groups is small, indicating convergence of the diagram- 
matic expansion. The CSB effect from this class is essentially negligible. 



The total CSB difference of the singlet scattering length caused by nucleon mass splitting 
amounts to 1.508 fm |]2D| which agrees well with the empirical value 1.6 ± 0.6 fm. Thus, 
nucleon mass splitting alone can explain the entire empirical CSB of the singlet scattering 
length. 

Starting from the CD-Bonn pp potential |2^, the parameters of the scalar-isoscalar 



bosons of that model have been adjusted such that the microscopically evaluated phase 
shift differences due to nucleon mass splitting are reproduced accurately. This yields the 
CD-Bonn nn potential. 



C. Phenomenological model 



An excellent example for a phenomenological construction of the CSB nuclear force is 
the recent Argonne Vig NN potential . As usual, the Argonne pp potential is fixed by a 
best-fit to the pp data. The Argonne nn potential is then constructed by starting from the 
pp one and readjusting the central force in the = 0, T = 1 state such that the empirical 
value for a^n is reproduced. For the 5* = 1,T = 1 state, where empirical information on nn 
scattering is not available, it is assumed that the CSB splitting of the central force is the 
same as in S* = 0, T = 1. 



D. Comparing the predictions for the two- nucleon system 

In Table |, we show the ^5*0 effective range parameters as calculated by the three models 
applied in this study. By construction, the p-u potential produces the same CSB difference 
as the nucleon mass (Mat) splitting model, namely, AacsB = 1-508 fm. The Argonne Vis 
potential yields AacsB = 1-654 fm. Thus, all three models have nearly identical results for 
AacsB which is exactly what we want as the starting point of our study. 

We now turn to energies Tiab > and calculate the CSB phase shift differences A6csb = 
Snn — Spp (without electromagnetic interactions) for all three models (see Fig. ^). In the ^5*0 
state at low energies, we have, of course, nearly identical phase shift differences because of 
the agreement in AacsB- However, as the energy increases, differences between the model 
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predictions emerge. The p-uo model maintains the largest A5csb above 150 MeV which may 
be explained by the fact that p-uo exchange is of shorter range than the 2tt exchange which 
the M^r splitting CSB potential is based upon. The differences in can be explained with 
the same argument. 

The largest differences between the model predictions occur in the ^Pj waves (cf. Fig. 
As expected, the p-uj potential now clearly reveals its large spin-orbit component typical for 
vector-meson exchange. Note that the spin-orbit force of ^Vp^ is of opposite sign as the one 
of ordinary one-omega exchange [cf. Eqs. ( p.9| ) and ( |2.3| ) and keep in mind that {p^\Hm\uj) 
is negative]. 

The Argonne Vi^ potential follows in ^Sq the trend of the p-uj mechanism and in ^Pq it 
is close to the MAr-splitting model. In the other partial waves, it is not close to any of the 
microscopic models for CSB. 

In summary, in spite of identical AacsB-, the p-uJ and the M^r-splitting models show 
drastic differences in ^Pj waves. Unfortunately, we do not have any empirical information 
on A6csB and, so, there is no direct way to tell which is right and which is wrong. Apart 
from the ^5*0 scattering lengths, the only empirical information on CSB that we have are 
some binding energy differences, to which P waves do contribute. Therefore, we will turn 
in the next section to such binding energy differences with the hope that the differences in 
P waves may impress a detectable signature. 



III. PREDICTIONS FOR SYSTEMS WITH A>2 

A. '^H-'^He binding energy difference 

The experimental value for the difference between the binding energies of and '^He is 
764 keV ('^H is more bound). Most of this difference is due to the static Coulomb energy 
(amongst finite-size protons) which accounts for 648 ± 4 keV |3^-0. Another 35 ± 3 keV 



come from electromagnetic effects neglected in the static Coulomb approximation 



and 14 ±2 keV are due to the n — p mass difference in the kinetic energy ||36|. After all these 
obvious corrections, a binding energy difference of 67 it 9 keV remains which is commonly 
attributed to CSB of the nuclear force. 

We have applied the three different CSB forces presented in the previous section in 



accurate momentum-space Faddeev calculations of the three-nucleon systems [32|. Our 
results for the ^H-^He binding energy differences are shown in the upper part of Table ||. 
We conducted two types of calculations. In the first type, we included CSB only in ^5*0 while 
all other partial waves are treated charge-symmetric. In the second type, CSB was included 
in all T = 1 partial waves (i. e., distinct pp and nn potentials were used in the T = 1 states). 

The predictions by the CSB models are between 60 and 66 keV. Thus, they are all 
consistent with the empirical value of 67±9 keV and no discrimination is possible. Moreover, 
the CSB contribution beyond ^5*0 is small, 2-5 keV (about 6% of the total), which is within 
the empirical uncertainty. Therefore, it is impossible to draw any conclusions concerning 
the CSB contributions from ^Pj waves. 

The trends in the results are consistent with the phase shift differences shown in Fig. ^: 
The p-u! model generates more binding energy difference from the ^5*0 state (3.3 keV more) 
and more from ^Pj waves (2.5 keV more) as compared to the M^r-splitting model. 
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B. Nolen-SchifFer anomaly 



It is a well-known experimental fact that the single-particle energies of corresponding 
states in mirror nuclei are different. If one assumes that the strong part of the nuclear force 
is charge symmetric, i.e. the strong proton-proton interaction is identical to the interaction 
between two neutrons, then these differences would originate entirely from the electromag- 
netic interaction (mainly Coulomb) between the nucleons. For this reason, it is customary 
to call these single-particle energy differences the Coulomb displacement energies. After 
accurate experimental data on the charge distribution became available from electron scat- 
tering experiments, Hartree-Fock calculations with phenomenological models for the NN 
interaction like the Skyrme forces were performed which reproduced these measured charge 
distributions with good accuracy. The Coulomb displacement energies which were evaluated 
with these Hartree-Fock wave functions, however, underestimated the experimental data by 
typically seven percent. This has become known as the Nolen-Schiffer (NS) anomaly ||2(J[| . 
Many attempts have been made to explain this discrepancy by the inclusion of electromag- 
netic corrections, many-body correlations beyond the Hartree-Fock approximation, or by 
explicit charge-symmetry breaking terms in the NN interaction P,|37|-^. In these inves- 
tigations it turned out that the CSB of the nuclear force (with the nn interaction more 
attractive than the pp one) is probably the major reason for the 'anomaly'. 

We will now study the impact of CSB of the A^A^ interaction on the single-particle 
energies for protons and neutrons in nuclear matter as well as in finite nuclei. 



1. Nuclear matter 

We calculate the single-particle potential for protons, Up{k), and neutrons, Un{k), 
as a function of the momentum k in symmetric nuclear matter using the self-consistent 
Brueckner-Hartree-Fock approach [|19|. From this we obtain the CSB energy difference 

AUcsB{k) = Up{k)-Un{k). (3.1) 

Since the momentum dependence of AUcsb is weak, we choose k = kp- Our results at 
nuclear matter density, kp = 1.35 fm~^, are shown in the lower part of Table f|. The 
most encouraging aspect of the results is that, here, we encounter a large contribution from 
states with L > (about 50% of the total). Consequently, we also observe a substantial 
difference between the model predictions, with the p-uj model producing about 20% more 
energy difference than the MAr-splitting model. 

Unfortunately, no reliable empirical estimates for AUcsb in nuclear matter exist such 
that we cannot draw any conclusions. Accurate data exist for finite nuclei which we will 
consider in the next section. 



2. Finite nuclei 

We choose ^^O for our sample nucleus and we wish to calculate Coulomb displacement 
energies around this nucleus HTJ. Calculations of this kind are very involved and, therefore. 
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we need to discuss first how to conduct such microscopic nuclear structure calculations in a 
proper way. 

One possibility would be to perform self-consistent Brueckner-Hartree-Fock (BHF) cal- 
culations and extract the Coulomb displacement energies from the single-particle energies 
for protons and neutrons. We do not take this approach, for the following reasons: (i) Such 
self-consistent BHF calculations typically predict the radii for the charge-density distribu- 
tions too small [^|. This implies that the leading Coulomb contribution to the displacement 



energy would be overestimated. Also the calculation of the correction terms would be based 
on single-particle wavefunctions which are localized too much, (ii) BHF calculations are 
appropriate for short-range correlations. However, long-range correlations involving the ad- 
mixture of configurations with low excitation energies in the uncorrelated shell-model basis 
require a more careful treatment, (iii) The BHF single-particle energies do not account for 
any distribution of the single-particle strength consistent with realistic spectral functions. 

For the reasons listed, we take the following approach. We use single-particle wave func- 
tions from Hartree-Fock calculations with effective nuclear forces, which yield a good fit to 
the empirical charge distribution. These wave functions are used to determine the leading 
Coulomb contribution and corrections like the effects of finite proton size, the electromag- 
netic spin-orbit interaction, the kinetic energy correction due to the mass difference between 
proton and neutron, and the effects of vacuum polarization. Actually, for these contributions 
we use the results by Sato . The first column of our Table |Tl| is taken from table 2 of 



Ref. [§8| which includes all the effects just mentioned. 

The correlation effects are taken into account in a two-step procedure. We assume a 
model space defined in terms of shell-model configurations including oscillator single-particle 
states up to the IpOf shell. We use the oscillator parameter b = 1.76 fm which is appropriate 
for ^^O. The effects of short-range correlations are calculated by employing an effective 
interaction, i.e. a ^-matrix suitable for the model space. This ^-matrix is determined as the 
solution of the Bethe-Goldstone equation 

= V + V ^ ^^Qd ^(^) ^ (3,2) 

" Vmod-' *^mod 

where T is identified with the kinetic energy operator, while V stands for the bare two- 
body interaction including the Coulomb interaction and accounting for CSB in the strong 



interaction. The Pauli operator Qj^qJ in this Bethe-Goldstone Eq. (3^) is defined in terms 
of two-particle harmonic oscillator states \a(3 > by 

{0 if a or /5 from Os or Op shell 
if a and (3 from IsOd or IpOf shell (3.3) 
> elsewhere 

As a first approximation we use the resulting ^-matrix elements and evaluate single-particle 
energies in the BHF approximation ea- This approximation, which will be denoted as BHF 
in the discussion below, accounts for short-range correlations, which are described in terms of 
configurations outside our model space. In a next step we add to this BHF definition of the 
nucleon self-energy the irreducible terms of second order in Q which account for intermediate 
two-particle one-hole and one-particle two-hole configurations within the model-space 
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^(2) _ 1 < cth\Q\piP2 >< piP2\G\ah > ^ 1 < ap\Q\hih2 >< hih2\G\ap > 



(3.4) 

Applying the techniques described in [^], we can solve the Dyson equation for the single- 
particle Greens function Gai^oj) 

Go.{u) = + go,{u)U'i\uj)G^{uj) (3.5) 

with Qa the BHF approximation for the single-particle Greens function, and determine its 
Lehmann representation 



2 



= E , z.^ , + E -^^^^ EmZT^— • (3-6) 



2 



This yields directly the energies of the states with A ± 1 nucleons that we are interested in. 

Our results for the Coulomb displacement energies are listed in Table |ITI| for various one- 
hole and one-particle states relative to ^^O. The first column of this table, C'^^\ contains the 
results of Ref. ||38| for the leading Coulomb contributions, the corrections due to the finite 



proton size, the electromagnetic spin-orbit interaction, the kinetic energy correction due to 
nucleon mass splitting, and the effects of vacuum polarization. As discussed above, we think 
that it is more realistic to evaluate these contributions for single-particle wave functions 
which are derived from Hartree-Fock calculations with phenomenological forces rather than 
using the wavefunctions derived from a microscopic BHF calculation. 

The second and third columns of Table ^ list the corrections to the Coulomb dis- 
placement energies which originate from the treatment of short-range {^sr) and long-range 
correlations {^lr) discussed above. The correction A^r has been derived from the dif- 
ferences of BHF single-particle energies for protons and neutrons subtracting the Coulomb 
displacement energy evaluated in the mean-field approximation 

^SR = ef^^ (proton) - ef^^ (neutron) - A^^^^^ ggj^ (3.7) 

In this case the BHF calculations have been performed with the np versions of the different 
interactions. The correction terms Am have been evaluated in a similar way from the 
quasiparticle energies determined in the Greens function approach, subtracting the BHF 
effects already contained in Asr. The correction terms Asr and A^r include the effects 
represented by irreducible diagrams of second and higher order in the interaction, in which at 
least one of the interaction lines represents the Coulomb interaction. In addition they contain 
the effects of folded diagrams discussed by Tam et al. [^. We find that the correlation effects 
are rather weak. The short- and long-range contributions tend to cancel each other. This is 
true in particular for the one- hole states Pg^g p^yg- The effects of short-range correlations 
dominate in the case of the particle states, d^/2 and lsi/2, leading to a total correlation effect 
in the order of 100 keV in the Coulomb displacement energies. This effect is slightly larger 
for the Argonne potential than for the Bonn potentials because of the stronger correlations 
in the case of Argonne. 
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The contributions to the Coulomb displacement energies caused by CSB of the A^A^ 
interactions, ^csB) are listed in the fourth column of Table III. We have conducted separate 
calculations for each of the three models for CSB introduced in Sec. II. 

We note that, also in the calculations of the Coulomb displacement energies, it is im- 
portant to include CSB beyond the ^5*0 state in the A^A^ interaction. Similarly to what we 
found for the nuclear matter AUcsb (cf. Sec. III.B.l and lower part of Table H), CSB in P 
and higher partial waves contributes about 50% of the total Acsb- This was demonstrated 
in Ref. ||. 

We achieve satisfactory or even good agreement in some cases, like p^^^ lsi/2, but 
there are discrepancies in others. For p^J^^ ^-nd ^5/2 the remaining discrepancies are larger 
than the Aqsb contribution. 

The general trend in the results is that the p-u model generates about 10-20% more 
Acsb than the MTv-splitting mechanism, and Argonne Vig is, in general, in-between the 
two. Even though the p-u trend is a favourable one, it is not sufficiently pronounced such 
that one could give a preference to this model. In the critical cases, like p'^^^ ^5/2, the 
discrepancies between all model predictions, on the one hand, and experiment, on the other, 
are much larger than the differences within the model predictions. 

Concerning the remaining discrepancy, a comment is in place. Note that the nuclear 
structure part of our calculations may carry some uncertainty. This is true in particular for 
the evaluation of the leading Coulomb contribution, which is sensitive to the Hartree-Fock 
wave functions which are used. To obtain an idea of how large such uncertainties may be, 
we compare the results for Coulomb displacement energies using the Skyrme II force and no 
CSB by Sato [3S] with the more recent ones by Suzuki et al. [P]. For the single-hole state 
Pi/^2^ Suzuki's result is larger by 167 keV as compared to Sato; and for the single-particle 
state d5/2, the two calculations differ by 138 keV. Uncertainties of this size can well explain 
the remaining discrepancies in our results. 



IV. SUMMARY AND CONCLUSIONS 

We have tested three different models for the charge symmetry breaking (CSB) of the 
A^A^ interaction. The models are based upon p-oo mixing, nucleon mass (Mat) splitting, and 
phenomenology (Argonne Vis). All models reproduce the empirical value for the CSB of the 
""^^o scattering length (AacsB) accurately. 

We reveal that there are considerable differences in the predictions by these models for 
CSB in ^Pj waves. We have investigated the impact of these differences on some observable 
quantities of A> 3 nuclear systems that are sensitive to CSB of the nuclear force. 

We find that the '^H-^He binding energy difference is essentially ruled by AacsB and that 
P and higher partial waves contribute only about 6%. Therefore, this quantity is unsuitable 
to discriminate between different models for CSB of the nuclear force. 

A test calculation conducted for nuclear matter shows that in heavier nuclear systems 
the difference between proton and neutron single-particle energies receives about 50% from 
partial waves other than ^5*0. Motivated by this result, we have calculated the Coulomb 
displacement energies around the closed shell nucleus ^^O. We find that the contribution to 
these energy differences from CSB of the A^A^ interaction (which is in the order of 100 keV) 
differs by 10-20% among the three models for CSB, which is appreciable. Unfortunately, 
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the nuclear structure part of these calculations, in particular the evaluation of the leading 
Coulomb contributions, carry an uncertainty in the order of 100 keV such that the subtle 
differences between the competing CSB models get swamped. Therefore we must conclude 
that, based upon the calculations conducted in this study, we are unable to give preference 
to any of the three CSB models. 

What we need are observables for which the nuclear structure part is fully under control. 
This suggests to look into nuclear few-body reactions for which exact calculations can be 
performed An example could be the analyzing power. Ay, in nucleon-deuteron scatter- 
ing that is known to depend sensitively on the ^Pj waves of the NN potential. Accurate data 
on p — d and n — d Ay exist and these data exhibit a clear signature of CSB. Unfortunately, 
these data cannot be explained without resource to three-body forces which may obscure 
the CSB aspect of the problem. In any case, we like to encourage the nuclear few-body 
community to identify spin observables in few-body reactions that show sensitivity to 
CSB of the NN interaction, particularly, the one that comes from ^Pj waves. Investigations 
of this kind may ultimately allow to discriminate between competing models for CSB of the 
nuclear force. 
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TABLES 



TABLE L ^5*0 scattering length (a) and effective range (r), for proton-proton (pp) and neu- 
tron-neutron (nn), with Coulomb effects (C) and without any electromagnetic effects (A^), in units 
of fm. 



p-uj mixing 



Mjv splitting 



Argonne Vi 



18 



Experiment 



,c 
pp 

1^ 
-^pp 

„N 
PP 

2^ 
nn 



^acsB 
^rcsB 



-7.8154 
2.773 

-17.460 
2.845 

-18.968 
2.816 
1.508 
0.029 



-7.8154 
2.773 

-17.460 
2.845 

-18.968 
2.819 
1.508 
0.026 



-7.8138 
2.787 

-17.164 
2.865 

-18.818 
2.834 
1.654 
0.031 



-7.8149 ± 0.0029'' 
2.769 ± 0.014'' 



-18.9 ±0.4'' 
2.75±0.11'= 

1.6 ±0.6 
0.10 ±0.12 



"Reference 
''Reference 
'^Reference fHl. 



TABLE II. Applications of CSB potentials in the three-nucleon system and in symmetric nu- 
clear matter. 









p-uj mixing splitting Arj 


^onne Vis 


Empirical 








'^H-'^He binding energy difference (keV) 






CSB 


in 


only 


60.9 57.6 


62.1 




CSB 


in 


all r = 1 states 


65.8 60.0 


65.1 


67 ± 9 








Nuclear matter AUcsB (MeV) 






CSB 


in 


only 


0.168 0.154 


0.180 




CSB 


in 


all T = 1 states 


0.367 0.311 


0.301 
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TABLE III. Coulomb displacement energies for single-hole ^'1/2) single-particle 

states ((^5/2 ^-ncl lsi/2) around ^^O. The single-particle contribution, C^^\ is from Sato |3^. Contri- 
butions due to short-range correlations, A^^, long-range correlations inside the model space. Am, 
and due to CSB of the strong interaction, AcsB, are calculated for three different CSB models. 
The total results for the displacement energies, are compared to the experimental data given 

in the last column. All entries are in keV. 





> 


A 

AsR 


A 

Alr 


i^CSB 




Exp 


P3/2 


p-to mixing 


3205 


-44 


46 


106 


3313 


3395 


Mjv splitting 




-44 


46 


97 


3303 






Argonne Vi^ 




-71 


47 


108 


3285 




Pl/2 


p-io mixing 


3235 


-52 


37 


124 


3344 


3542 


Mjv splitting 




-52 


37 


102 


3322 






Argonne Vig 




-79 


39 


103 


3297 




^5/2 


p-oj mixing 


3135 


154 


-15 


93 


3367 


3542 


Mtv splitting 




154 


-15 


87 


3361 






Argonne Vig 




187 


-18 


92 


3401 




lSl/2 


p-to mixing 


2905 


159 


-45 


134 


3154 


3166 


Mtv splitting 




160 


-46 


112 


3132 






Argonne Vig 




198 


-47 


112 


3174 
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FIGURES 



FIG. 1. p°-a; exchange contributions to the nuclear force that violate charge symmetry. 



-it", a' , p^, to, S 



P P 

(a) 



(b) 



FIG. 2. One-boson-exchange (OBE) contributions to (a) pp and (b) nn scattering 



19 



n A+ 



P 



A 



P A 



(a) 



A" 



n A- 



V 



p A« 



V 



n A-" 



(b) 



FIG. 3. Two-pion-exchange contributions with NA intermediate states to (a) pp and (b) nn 
scattering. 
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FIG. 4. CSB phase shift differences (5„„ — Spp (without electromagnetic interactions) for lab- 
oratory kinetic energies Tiab below 350 MeV and partial waves with L < 2 as generated by p-uj 
mixing (solid line) , nucleon mass splitting (dashed) , and the phenomenological Argonne Vis model 
(dotted). 
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